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ABSTRACT 

Using Hubble Space Telescope (HST) photometry, we age-date 59 supernova remnants (SNRs) in 
the spiral galaxy M31 and use these ages to estimate zero-age main sequence masses (Mzams) for 
their progenitors. To accomplish this, we create color-magnitude diagrams (CMDs) and employ CMD 
fitting to measure the recent star formation history (SFH) of the regions surrounding cataloged SNR 
sites. We identify any young coeval population that likely produced the progenitor star, then assign 
an age and uncertainty to that population. Application of stellar evolution models allows us to infer 
the MzAMsfrom this age. Because our technique is not contingent on identification or precise location 
of the progenitor star, it can be applied to the location of any known SNR. We identify significant 
young star formation around 53 of the 59 SNRs and assign progenitor masses to these, representing a 
factor of ^2 increase over currently measured progenitor masses. We consider the remaining 6 SNRs 
as either probable Type la candidates or the result of core-collapse progenitors that have escaped 
their birth sites. In general, the distribution of recovered progenitor masses is bottom heavy, showing 
a paucity of the most massive stars. If we assume a single power law distribution, dN/dM oc M a , 
we find a distribution that is steeper than a Salpeter IMF (a = —2.35). In particular, we find values 
of a outside the range —2.7 > a > —4.4 to be inconsistent with our measured distribution at 95% 
confidence. If instead we assume a distribution that follows a Salpeter IMF up to some maximum 
mass, we find that values of Mmux > 26 are inconsistent with the measured distribution at 95% 
confidence. In cither scenario, the data suggest that some fraction of massive stars may not explode. 
The result is preliminary and requires more SNRs and further analysis. In addition, we use our 
distribution to estimate a minimum mass for core collapse between 7.0 and 7.8 M Q . 
Subject headings: galaxies:individual:(M31) — supernovae: general 



1. INTRODUCTION 

Both theoretically and observationally, core-collapse 
supernova (CCSN) explosions are linked with the death 
of massive stars. However, because precise measurements 
of CCSN progenitor masses are scarce, the mapping 
between explosion scenario and the progenitor mass 
distribution is less clear. The most common method 
of progenitor mass determination is direct imaging, in 
which one identifies the progenitor star in multi-band 
pre-explosion imaging, fits a spectral energy distribution 
to the photometry, and verifies the star is gone once 
the SN has faded. This methodology is ideal given 
appropriate data, but it is limited to contemporary SNe 
that have pre-explosion Hubble Space Telescope (HST) 
or deep ground-based imaging. As a result, only ^25 
SNe have any constraint on their pro genitor masses, and 
half of these are only upper limits ([Smartt et all 120021: 
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in the number of progenitor masses could allow us a 
powerful window into further understanding supernovae. 

Direct imaging has created a prototypical picture of 
massive star death. Type II-P SNe are assumed to be 
created by red supergiant (RSG) stars with intact hy- 
drogen envelopes, while more exotic SNe may be cre- 
ated by higher mass stars. However, a number of ques- 
tions remain to be a nswered regarding supernova physics. 
iSmartt et all (|2009f ) identify what they term to be the 
'red supergiant problem,' an observed lack of progenitors 
of Mzams~16-30 M Q that we would expect to explode as 
Type II-P SN. A variety of channels have been proposed 
to account for these missing explosions, including direct 
black hole formation or the stars exploding as different 
SN types. Binarity of the progenitor system can also play 
a significant role in SN explosions. iSmith et all (|2011al ) 
show that numbers of Type Ibc and lib explosions are 
significantly underestimated if single-star explosions are 
assumed to be the only SN channel. 

These issues arc difficult to address without further 
progenitors to analyze. Unfortunately, while direct- 
imaging of the progenitor star would be ideal, there are 
constraints that limit the frequency with which it may 
be applied. First, direct imaging is limited by the SN 
rate: roughly one per century in a large spiral galaxy 
(jCappellaro et al.lll999l ). Second, direct imaging requires 
that the host galaxy be close enough to r esolve the pro - 
genitor star, an effective limit of ^20 Mpc (Smartt 2009). 
Third, preexisting sub-arcsecond images of the site must 
exist, necessitating archival HST images with either the 
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Wide Field Planetary Camera 2 (WFPC-2) or the Ad- 
vanced Camera for Surveys (ACS). Finally, even with 
the existence of appropriate images, there is no guar- 
ante e that the pr e curso r will be identified. For exam- 
ple, iSmartt et al.1 (|2009( ) summarized 20 progenitor de- 
tections in Type H-P explosions. Only five progenitors 
were directly observed, and an additional two which fell 
on compact star clusters were well constrained. The re- 
mainder had no detection, but were given upper bounds 
based on the maximum luminos ity with whi c h the y 
could escape detectio n (see also | Smartt et al.l (|2003f ): 
iVan Dvk et all IpOOSch : iMaund fc Smarttl (|2005m . 

In short, while direct imaging is ideal given appropriate 
data, it suffers from prerequisites that limit the number 
of opportunities in which it may be applied. It is clear 
that development of an independent and complementary 
technique would be of interest. 

In this paper we make use of an alternative technique, 
stellar population analysis, to infer progenitor masses. 
The technique involves examining a color-magnitude di- 
agram (CMD) of the surrounding population of stars to 
infer an age and mass for the progenitor star. The ad- 
vantage of our methodology is that we are not reliant 
on individual identification and photometry of the spe- 
cific progenitor star, leaving us free to apply our method 
to cataloged supernova remnants (SNRs) in addition to 
directly observed SNe. This ability drastically increases 
the number of progenitor masses we may find, allowing 
us to make a more complete measurement of the underly- 
ing distribution of progenitor masses. Our method offers 
no direct way to probe the binarity of the system, which 
is relevant given the role binarity likely plays in certain 
CCSN scenarios. However, our estimate of the age of the 
surrounding stellar population is not affected by the par- 
ticular details of a given progenitor system. Our method 
also offers no way to determine the type of CCSN explo- 
sion we are observing, so for SNR analysis we are not 
able to provide a direct link between SN explosion type 
and progenitor mass. However, we will be able to com- 
ment on the ranges of ZAMS masses that are producing 
CCSN explosions. 

Stellar population analysis has been used many times 
historically to analyze the characteristics of progenitors. 
Many groups have carried out age dating of stellar clus- 
ters coincident wit h supernoyae and derived correspond- 
ing stellar masses (lEfremovl 1 1 9911 : iWalborn et all If 9931: 
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(12009ft examined the SFH m ap of the LMC published 
by lHarris fc Zaritskvl (|2009i). They estimate ages and 
masses for the f our most recent SNR pro duced by CCSNe 
in the LMC. In iGogarten et all (|2009aft . we employed a 
technique identical to the one used in this paper to es- 
timate the age of the progenitor of NGC 300 OT2008-1. 
We found that the star was likely to have formed 8-13 
Myr ag o, corresponding to a Mzams of 12-17 M Q . Fi- 
nally, in lMurphv et al.l ( 20111 ) we applied the same tech- 
nique as this paper to archival images of SN 201 ldh, 
finding a most likely age of 17^4 Myr and a Mzams of 
13j^ Mq. Clearly there is a long list of successful applica- 
tions of stellar population analysis leading to constraints 
on progenitor physical properties. 



In this paper we examine 59 SNRs in the spiral galaxy 
M31, listed in Table [TJ In §2 we outline the methods 
used to determine star formation histories (SFHs), ages, 
and masses for the progenitor stars. In §3 we present 
the results of the analysis. We find that 53 SNR display 
recent star formation, and we derive masses for these. 
We also present four example cases that are representa- 
tive of the overall sample. Finally, in §4 we examine the 
distribution of progenitor masses. We search for indica- 
tions of a minimum mass, which we find to be between 
7.0 and 7.8 Mq, and note that the distribution is more 
bottom-heavy than a Salpcter IMF distribution. 

2. METHODOLOGY 

The basic procedure for deriving a mass estimate for 
each progenitor is as follows. We first use coordinates 
from three SNR catalogs as explained in §2.2, which 
we assume are accurate to within a few arcseconds and 
unbiased to any one type of SN. We then search for 
any appropriate HST fields that contain the SNR in 
question. We perform photometry on the stars in each 
field, and create CMDs of the region within ~50 pc of 
the SNR coordina tes, assuming a distance m odulus of 
m — M = 24.47 (jStanek fc Garnavichl [1998ft . We use 
CMD fitting to measure the SFH of each region. For 
each region that displays recent SF, we use the age of this 
recent SF to associate an age with the progenitor star. 
Finally, we apply stellar evolution models to convert the 
age of the progenitor to the Mzams of the progenitor. 

2.1. Key Assumptions 

Our method requires several assumptions to produce 
results. Above all, we must assume that some of the 
stars surrounding the SNR are in fact coeval with the 
progenitor star. There are several pieces of evidence 
suggesting this assumption is reasonable. First, 90% of 
stars form in clusters with sizes of ~1 pc (|Lada &: Ladal 
12003ft . These stars are expected to stay spatially asso- 
ciated for timcscalcs greater than the lifetimes of CC- 
SNe progenitors, even for dissolving, unbound clusters. 
Second, theoretical predictions are that stars will stay 
spatially associated on sca les of —100 pc for 100 Myr 
(|Bastian fc Goodwinl 120061) . Observ ational constraints 
su pport this (IGogarten et all l2009bT) . as do simulations 
bv lEldridee et al.l (|2011D ^ho found that -85% of CC- 
SNe will explode within 100 pc of their birth site. In 
addition, even if a progenitor would travel farther, we 
would still expect to see a fraction of the coeval young 
stellar population which is sufficient for age dating. The 
precise number of interest is the fraction of young stars 
that will stay within a certain distance in the lifetime of 
a core collapse supernova progenitor. The question is a 
complex one, and we know of no research that quantifies 
it precisely at this time. We conclude that it is a rea- 
sonable assumption that most of the young stellar pop- 
ulation around our observed SNRs is c oeval with their 
progen itors. Based on the analysis in IGogarten et al.l 
(|2009aft . we adopt a value of ~50 pc for the radius of our 
star selection annulus. 

Our method docs not assume any information about 
the type of SN that exploded. We assume that SNRs 
associated with recent star formation are CCSNe, but 
do not distinguish individual subtypes. For CCSNe, the 
particular type of supernova has no effect on our mass 
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determination process. However, the progenitors of ther- 
monuclear Type la SNe are likely to arise from older 
stellar populations which may vary considerably more in 
age than the coeval populations of massive stars. While 
there exists considerably discussion in the literature as 
to the precise nature of Type la progenitor systems and 
their distribution in age (see iMaoz fe Mannuccil (|2012f ) 
and references therein), we would expect to see Type la 
progenitors ranging from ages of a few hundred Myr to 
several Gyr. Since these timescales are comparable to 
or significantly greater than the dynamical timescale of 
M31, our methodology of identifying and analyzing co- 
eval stars is ineffective at measuring the precise ages of 
Type la progenitors. 

Our method does have the risk of misinterpreting a 
Type la SNR as a CCSNR. I n the ir volume limited sam- 
ple (D < 60 Mpc), iLi et al.1 (poTLT ) find that Type la SN 
compose 24% of observed SN. However, there are several 
qualitative reasons as to why one would expect a smaller 
Ia/CCSN fraction in our survey. 

First, there is some discussion as to whether such a 
volume-limited survey may underestimate the fraction of 
CCSN due to t heir faintness compared t o Type la SN (see 
discus si ons in iThompson et al.l (|2009[ ): iHoriuchi et al.l 
(|2009t) ; IHoriuchi et al.l (|201lD ). This potential bias sug- 
gests that 24% could be considered an upper limit. In 
addition, the areas of M31 with greater HST coverage are 
primarily star forming regions, where we would naively 
expect CCSN to be more common, increasing the rela- 
tive fraction of CCSN. Thus the actual fraction of Type 
la SNR we expect in our sample is probably somewhat 
less than 24%. 

Moreover, it is likely that some Type la SNR can be 
identifie d by their lack o f an as sociated young stellar pop- 
ulation. iBadenes et al.1 (|2009[ ) examined four la sites in 
the LMC and found that three displayed a significant 
lack of recent star formation. Out of the 59 examined, we 
found 6 SNR displayed no recent SF in their surrounding 
stellar populations. 

Finally, we have examined what effect additional hid- 
den la contamination will have on our data. The proce- 
dure we used is explained in §4.3. We found that addi- 
tional la contamination up to a total fraction of ~25% 
has very little effect on the overall distribution observed 
in the sample. 

To summarize, we argue that there are several quali- 
tative reasons to assume that the effective la fraction in 
our sample is smaller than 24%. In addition, our method 
offers a means to identify and remove some SNR that 
arc likely la progenitors. Finally, our overall results are 
largely insensitive to inclusion of additional la contam- 
ination up to the 24% obs erved in the volume-limited 
survey from ILi et all (pOll . 

Another possibility for SNRs with no coincident young 
SF is that they are the result of high- velocity progeni- 
tor stars that have left their birth sites. As a test of this 
possibility, we examined any site with no observed recent 
SF in GALEX FUV data. We found that all 6 sites with 
no recent SF are still relatively close (~50 to 100 pc) to 
GALEX FUV sources. If we assume that the progeni- 
tor star resulted from the most nearby FUV source for 
each zero-SF SNR, then a 50 Myr old star would only re- 
quire velocities of a few km/s to reach the SNR site from 
the FUV source site, well within reasonable velocities for 



runaway stars. Higher mass stars would of course require 
higher velocities to leave the FUV source in their shorter 
lifetimes. This analysis suggests that our study is unable 
to differentiate between Type la SNR and SNR resulting 
from high-velocity stars. We note, however, that a null- 
result for progenitor mass will not affect the distribution 
of progenitor masses we measure. It will simply reduce 
our effective sample size. 

Finally, our method is highly contingent on the accu- 
racy of the stellar evolution models used to model our ob- 
served CMDs. Because the models are theoretical, they 
generally do not have easily quantifiable uncertainties. 
The wide array of available model sets tend to system- 
atically differ in bolomctric luminosity and temperature, 
but are consistent within about 0.2 mag in bolometric 
luminosity and 0.02 dex in log temperature. We there- 
fore quantify the uncertainty of our results due to the 
models used by including random shifts in the bolomet- 
ric luminosity and temperature of the models a s part of 
our Monte Carlo (MC) tests (see also iDolphir] (|2002D ). 
We detail the process further in §2.5. 

2.2. Data and Photometry 

There are only a few ext ensive catalogs of SNRs i n 
M31. We use the cata logs ofjB raun fe Walterbosl (fl99lh . 
iMagnier et all (fl995l ). and iWilliams et al.l (|1995f ). All 
thre e make identification s based on [SII]-to-Ha ratios, 
and IMagnier et al.1 (|1995| ) make a dditional use of mor - 
phology and OB star associations. IMagnier et al.1 (|1995l ) 
identify three confidence levels based on approximate lev- 
els of contamination, although they note that the con- 
fidence levels are somewhat subjective. We only per- 
form analysis on their first and second confidence can- 
didates. The catalogs are not mutually exclusive, and 
the coverage areas of the three catalogs overlap in sev- 
eral areas, causing several remnants to be identified in 
both catalogs. To avoid treating a double-identified 
remnant as two separate remnants, we identify all rem- 
nants from two separate catalogs in which the 50 pc ar- 
eas overlap and only select coordinates from one cata- 
log. In these d uplicate cases, we a d opted the coordi- 
nates first from iBraun fe Walterbosl (119931) . then from 
IMagnier et all (|1995l ) if IB raun fe Walterbosl (fl993l) did 
not identify a remnant in that position. On occasion, a 
single catalog identified two remnants in which these 50 
pc regions overlapped, but we still performed analysis on 
both candidates in such a case. In addition, we cross- 
referenced the SNR catalogs with iSasaki et al.1 ()2012l ) 
and eliminated two, K567 and K884, which were found 
to not be SNR based on their optical/X-ray properties. 

For our analysis, we require that the stellar population 
surrounding the SNR be imaged by either the Advanced 
Camera for Surveys (ACS) or the Wide-Field Planetary 
Camera 2 (WFPC2) instruments on HST in at least two 
broadband filters. We also require our 50% complete- 
ness limits (see §2.3) for a given CMD to be at least 
24.5 magnitude in F475W, F555W, or F606W, whichever 
represents the blue filter (all qualifying fields of suffi- 
cient depth were either F475W, F555W, or F606W vs. 
F814W). The location of the MS turnoff is very similar in 
color in these filter sets, so adopting a similar magnitude 
cut in all three is reasonable. 

Our comparisons between SNRs imaged multiple times 
at varying depths indicate that simply reaching the MS 



4 



turnoff for a given CMD is insufficient to correctly deter- 
mine the age. CMD fitting is highly sensitive to densities 
of stars at given ages, requiring a well-sampled stellar 
population at the age of interest. In general, we found 
that our methodology applied to shallower data had a 
tendency to miss older bursts of SF which still may have 
resulted in CCSN. The available 2-filter broadband data 
for M31 happens to be distributed in two groups, includ- 
ing many shallow fields (blue depth < 23) and many deep 
fields (blue depth > 25), with very few in between. The 
choice of our depth cut causes us to reject five SNR that 
have only been imaged with shallow WFPC-2 data. 

In Table[T]we list the SNR for which our criteria is met, 
along with the corresponding HST fields. In cases where 
the regions in question were imaged multiple times, we 
selected whichever set of images had the greatest number 
of stellar detections. In Fig. [1] we chart the locations 
of the examined remnants on a star-subtracted image 
of M31 in Ha. Those remnants colored red represent 
probable la candidates, where no young SF was found. 
There is a clear grouping of CCSN candidates along the 
star forming arms of M31. 

We performed resolved stellar photometry using the 
photometry pipeline deve loped for the ACS N earby 
Galaxy Treasury program ([Dalcanton et al.l 120091 ). This 
pipeline uses th e DOLPHOT stellar photometry package 
(|Dolphinl 120001 ) to fit the well-characterized ACS point 
spread function to all of the point sources in the images. 
We then converted fluxes to Vega magnitudes using the 
standard zero-points and aperture corrections from the 
ACS handbook. We assess photometric errors and com- 
pleteness using fake star tests. At least 10 5 tests are 
performed by inserting fake stars of known color and 
magnitude into the data one at a time and blindly at- 
tempting to recover them with the same software. Both 
fake star tests and photometry were performed on the 
full HST fields. 

2.3. Calculation of Star Formation Histories 

CMD fitting is a powerful tool for measuring star for- 
mation histories (iGallart et al J 120051: iBarker et al.ll2007i : 
I Williams etHl l2009al lbl: iTolstov et al.ll2009f ). To esti- 
mate star formation histories of the regions surround- 
ing the SNR we used the software package MATCH 
()Dolphinll2002l) . MATCH works by creating many model 
CMDs based on theoretical isochrones for a variety of 
ages and metallicities. A linear combination of these 
model CMDs ar e then fit to t he ob serv ed CMD. We 
use th e models of iMarigo et al.l (|2008l ) and lGirardi et al.l 
(|2010ft both for CMD fitting and estimation of M ZA ms- 

For each field, we selected stars from our photome- 
try catalog within a 15" radius around the SNR co- 
ordinates, which equates to a physical size of ^50 
pc at an assumed distance mo dulus of m — M ~ 
24.47 (IStanek fc Garnavichl [l998h . This is consistent 
with the spatial correlation discussed in §2.1. We ac- 
count for photometric errors through the use of fake star 
tests. We selected fake stars in a region ~2.5 times the 
radius of the real star annulus to ensure that at least a 
few thousand recovered fake stars are included for each 
SNR. We defined lower magnitude limits as the point at 
which fake star completeness dipped below 50%. 

MATCH requires a variety of parameters to generate 
and fit CMDs. We assumed a Salpetcr IMF, dN/dM = 



M~ 2 35 (lSalpeterlll955|) . and a binary fraction of 0.35. 
iGogarten et al.l (|2009al) demonstrated that varying the 



IMF value from -2.0 to -2.7 or varying the binary frac- 
tion from 0.2 to 0.5 had no effect on the epoch assigned 
to a recent burst. We also varied the IMF from -1.3 to 
-3.3 for select regions and found no significant difference 
in relative fractions of SF for various epochs. Although 
absolute amplitudes of SF did change, only relative am- 
plitudes have an effect on our method of determining 
ages. Note that the IMF and binary fraction are used by 
MATCH purely for purposes of populating the models. 
We assume nothing about either value with regards to a 
potential CCSN progenitor system or the overall progen- 
itor mass distribution. 

MATCH produces fits in logarithmic age bins. We 
adopted 71 age bins increasing in .05 increments from 
6.60 (4 Myr) up to 10.10 (12.5 Gyr). We are unable to 
fit for ages younger than this due to the lack of isochrones 
at younger ages. MATCH will interpret any SF from a 
population younger than the 6.60 to 6.65 (4 to 4.5 Myr) 
bin as being included in this bin. Thus our youngest bin 
actually includes all SF from present times back to 4.5 
Myr. As a result, we may only quote an upper age limit 
for any SF found in this youngest bin. In Fig. [2] we plot 
the isochrones for every-other age bin from 6.60 through 
8.00 in F555W vs. F555W-F814W. 

Fitting mctallicity with MATCH is not viable given 
the very weak dependence of the optical colors of the 
upper main sequence on metallicity. We therefore con- 
strained the metallicity to a spread of ~0.15 dex and 
to increase with the lifetime of the galaxy, which con- 
sistently produced best-fit metallicities of solar. This 
value is consistent with the known gas-phase metal licity 
of M31 (jDennefeld fe Kunthlll98l : iBlair et al.lll982l ). Fi- 
nally, we binned our CMDs in units of 0.3 in magnitude 
and 0.15 in color. This binning accounted for the fact 
that some fields had few upper MS stars, allowing us to 
reduce the impact of this paucity on our fits. However, 
we still maintain a number of bins significantly larger 
than the number of free parameters used in the fitting. 
We experimented with using finer binning and found no 
change in our age determinations within our uncertain- 
ties. 

2.4. Treatment of Reddening 

Recent SFHs are very sensitive to the treatment of red- 
dening. Our most simplistic treatment of reddening is 
to assume that all reddening is due to a galactic fore- 
ground, search over a specified range of reddening val- 
ues for a best fit, and apply this value to the CMD as 
a whole. However, we frequently found significant dif- 
ferential reddening (dA v ) across our small SNR-centered 
regions, which manifests as an increased width to the MS 
and red clump. Star forming regions will typically have 
extensive amounts of dust and gas, often not uniformly 
distributed, and as such there is no reason to believe that 
a single reddening value would correctly describe the re- 
gion as a whole. MATCH by default allows for 0.5 mag of 
differential reddening for populations below 40 Myr, af- 
ter which it falls linearly down to 0.0 mag at 100 Myr. In 
addition, the user has the ability to add additional full- 
field differential reddening across the CMD. We almost 
always found that the default treatment of reddening was 
inadequate to account for all the reddening present in the 
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regions. 

As an example of how ignoring differential reddening 
can lead to an erroneous result, we consider the SNR 
BW-102. The photometry is quite deep, extending to 
magnitudes of F475W = 27.6 and F8UW = 26.8. In 
the left panel of Fig. |3] we plot the observed CMD of 
the region, as well as the best-fit model generated by 
MATCH in the background in greyscale. In the right 
panel, we plot the cumulative SFH as a fraction of to- 
tal SFH in the past 50 Myr, assuming a fixed differ- 
ential reddening of dA v = 0.5 mags affecting only the 
young stars. We found a best-fit reddening value of 
A v = 0.975. All the star formation is concentrated in 
the youngest time bin (< 4.5 Myr), corresponding to a 
lower limit Mzams of 52 M Q . However, several quali- 
ties of the CMD would cause us to be suspicious of this 
result. The main sequence appears very dim, and ap- 
pears to contain lack the young, bright, blue stars that 
the model predicts at F475W < 22.5. We would not 
expect this given a result of very extensive young star 
formation. Looking beyond the main sequence, we see 
that the red clump is extended over nearly 2 magnitudes 
in F475W. Indeed, the model very poorly models the 
stars at F475W - 25.5, F475W - F8UW ~ 1.75Thus 
not only was the default value of differential reddening 
insufficient to model the entire CMD, the treatment of 
differential reddening as something unique to the young 
stellar population is clearly not reflected in the reality 
of the CMD. For regions such as this, it is necessary to 
apply differential reddening across the CMD. We must 
then define the amount of differential reddening to add. 

To understand our technique for determining the 
amount of dA v to use, we must expand on the differen- 
tial reddening model used by MATCH. MATCH applies 
differential reddening to the model CMDs by applying a 
top-hat distribution along the reddening line. The low 
end of the top-hat is defined by a foreground value of A v 
for which MATCH fits. The width of the distribution is 
then defined by the specified value provided by the user. 
When applying additional differential reddening, fit val- 
ues tended to improve until we reached unphysical values 
of dA v The only true constraint we may apply is that this 
minimum value of A v must be at least equivalent to the 
Milky Way foreground value. iSchlegel et al.l ([1998) find a 
foreground reddening value to M31 of A v = 0.19. While 
this value is likely to change slightly over the angular size 
of M31, the changes will clearly be small compared to 
the magnitude of differential reddening internal to M31. 
Our solution to CMDs such as BW-102 is to increase the 
width of the differential reddening distribution applied 
to all ages until the best-fit value of A v MATCH finds 
reaches the foreground value. We then use this differ- 
ential reddening distribution to measure the SFH and 
derive the age of the SN progenitor. 

In Fig. 3] we plot the SFH results for a variety of differ- 
ential reddening values, increasing from dA v = 0.0 up to 
dA v = 1.5 magnitudes, in increments of 0.1. Note that 
this differential reddening is in addition to the default 
young dA v = 0.5, leading to up to 2.0 magnitudes of 
differential reddening on the very youngest stellar pop- 
ulations. We found that a value of dA v = 1.5 was the 
point at which the overall reddening value dropped to the 
expected foreground value for BW-102. For this value of 
differential reddening, the SFH corresponded to an age 



of 36 Myr and a mass of ^9 M©, a much more rea- 
sonable result given the observed CMD. Note that the 
SFH reaches this result (within uncertainty) before the 
best-fit A v value reaches foreground, suggesting that our 
result is not highly contingent on the precise value of fore- 
ground reddening assumed. In the left panel of Fig [5] we 
plot the observed CMD with the best-fit model plotted 
in greyscale in the background, created using the above 
dA v procedure. In the right panel, we plot the cumula- 
tive SFH of the region for the past 50 Myr, with error 
bars from the Monte Carlo analysis described in §2.5 

We found essentially all fields suffered from differen- 
tial reddening in excess of the default values used by 
MATCH. When no extra differential reddening was in- 
cluded in the models, many regions without obvious 
young star formation were best fit by high overall redden- 
ing values and included star formation in the youngest 
few time bin. The issue arises from faint main sequence 
stars revealed by deep photometry. Without inclusion 
of differential reddening, these stars are fit by the uni- 
form foreground reddening value for the region, which is 
typically a high value (e.g. A v > 0.6). Using this erro- 
neously high value, the stars are corrected to be brighter 
and bluer than they actually are, leading MATCH to 
mistake them as evidence of young star formation. By al- 
lowing for differential reddening, we can incorporate the 
full range of different reddening values across the field. 

For our final procedure, we increase the amount of full- 
CMD differential reddening until we find the best fit A v 
reaches the foreground value, Ay = 0.19. We then adopt 
this value of differential reddening for subsequent analy- 
sis of that particular SNR. 

2.5. Assessing Uncertainties 
2.5.1. Multiple Coeval Populations 

Our ultimate goal is to estimate the mass of the pro- 
genitor. We do this by identifying the burst of recent 
star formation from which the progenitor star most likely 
originated, allowing us to assign an age to each progen- 
itor star. Ideally, each SFH would have a single isolated 
burst, making it easy to associate a single age to the 
young stellar population. In many cases, however, mul- 
tiple bursts of varying rates and durations appear in the 
young star formation history. As a result, we examine 
the cumulative star formation history of each CMD as a 
fraction of the total recent star formation. When exam- 
ined as such, the fraction of star formation in each age 
bin corresponds to the probability that the progenitor 
will be of that age. Because we are only interested in 
recent SF, we examine only the most recent 50 Myr of 
SF (see §4.1 for an explanation of this limit). 

2.5.2. Fitting Errors 

Our results suffer from both random uncertainties due 
to the sampling of the CMD and systematic errors due 
to differences between the theoretical isochrones and ob- 
servations. For example, if the models were consistently 
redder than the observed stars, we would generally find 
younger ages when these models are applied to real data. 
The random errors are generally highly dependent on the 
number of upper main sequence stars in the field, with 
higher numbers providing a tighter constraint on the data 
(although more upper main sequence stars may also sim- 
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ply ind icate a younger population). See iGogarten et all 
i|2009al) for further discussion of this point. 

To analyze our uncertainties, we performed a series 
of Mont e Carl o (MC) realizations of the data (see also 
iDolphinKpOOl i lWeisz etlUlpoTl ). In each MC run, we 
re-sampled the observed CMD to account for the Poisson 
errors of the stars. In addition, to estimate how system- 
atic model differences can impact the results, we add in 
random shifts to the models for each run. We used ran- 
dom shifts of a = 0.02 in temperature and a = 0.17 in 
bolometric luminosity. These systematic offsets are much 
larger than uncertainty in the distance modulus, thereby 
incorporating both effects into our overall uncertainty 
analysis. 

Uncertainty in the age of the recent star formation 
manifests itself in two distinct ways. The first is uncer- 
tainty measured by the MC analysis performed above, 
while the second is due to the width of the intrinsic 
spread of the burst across multiple age bins. We esti- 
mate the latter of these as the difference between the 
median age and the ages where 16% and 84% of star for- 
mation has occurred. We estimate the former of these as 
the RMS difference between the median age of the best 
fit and the median ages of SF from the MC tests. We 
then add these differences in quadrature to assign a confi- 
dence interval to each result. Thus for each progenitor we 
determine a median age, as well as (potentially asymmet- 
ric) uncertainties about this median. Finally, we make 
the assumption that we are unable to determine the age 
of SF to a greater precision than the age bins in which we 
have measured this SF. We thus round the age range to 
the age of the next isochrone out from the median. We 
perform this step for both younger and older star forma- 
tion, always rounding away from the median. The ages 
of these isochrones determine our final uncertainties. 

In order to use this method of analysis, we must define 
the maximum age of star formation which may produce 
core-collapse progenitors. As explained in §4.1, we use 
our distribution to estimate a minimum mass for core- 
collapse between 7.8 and 7.0 M . As a result, we adopt 
50 Myr as the maximum age of a core-collapse progenitor 
and only perform the above mass estimation analysis over 
the most recent 50 Myr of SF. 

2.6. Converting SFH to Progenitor Mass 

The IMarigo et all (p00l and iGirardi et all (pOlOh 
models specify a maximum mass for an isochrone at a 
given age. More massive stars will have already died 
off. Thus the isochrones wc have identified as bounding 
our confidence interval can be linked directly to values 
for Mzams ■ Our value for the median progenitor mass 
comes from interpolating the final isochrone masses be- 
tween isochrones to the median age value. This is neces- 
sary because our median age won't line up exactly with 
a defined isochrone. 

In Fig. [51 we p lot the IMarigo et all ([20081 ) and 
IGirardi et al.l (|2010f) isochrones for final isochrone mass 
vs. age for metallicities of Z = 0.004, 0.008, 0.019, and 
0.030. The models produce very similar age to mass 
conversions regardless of the assumed metallicity. The 
vast majority of our regions produce best fit metallici- 
ties of approximately solar. As a result, we adopt the 
Z = 0.019 isochrone for mass determinations. Note 
that masses change very quickly for younger populations, 



while masses for older populations change at a much 
slower rate. This means that our results for less mas- 
sive progenitors will naturally be more precise. 

Note that wc have neglected systematic uncertain- 
ties in the age-to-mass conversion process for our in- 
dividual progenitor results. This leads to very small 
error bars on some progenitor masses, especially those 
at older ages, where mass doesn't change significantly 
as a function of age. Our mass distributions are 
not sensitive to this model-dependent systematic un- 
certainty, as shown by similarity between the median 
progenitor distribution and the probability distribu- 
tion (sec §4.2). However, uncertainties are almost 
certainly underestimated for some progenitors, espe- 
cially at the low-mass end. To estimate the mag- 
nitude of this systematic uncertain ty, we compared 
the s olar m etallicity isoc h rones of IPietrinferni et al.1 
(p00l to the IMarigo et all (poM ): IGirardi et al.1 (120101 ) 
isochrones used for our anal ysis. We found t h at for 
a -50 Myr lifetime star, the IPietrinferni et al.l (|2004f ) 
isochrones find a maximum mass of 6.7 M^, compare d 
to the^ 7.3 M^ pred icted for the IMarigo et all (|2008T ): 
IGirardi et al.l (120101 ) isochrones. At a lifetime of ~22 
Myr, the youngest stars tested by IPietrinferni et al.l 
(p00l . the' maximum masses for the different isochrone 
sets are 9.9 M Q and 10.9 M Q respectively. These suggest 
systematic uncertainties of around 0.5-1.0 M Q for the age- 
mass conversion process. Again, we do not include these 
systematic uncertainties in the reported values in Table 
2; the reader should keep this in mind when interpreting 
the results of any individual progenitor star in our study. 

3. RESULTS 

In Table [2] we list results for all 59 SNR analyzed. Wc 
tabulate the designation, mass, age, number of main 
sequence stars (defined as FA75W - F814PF < 0.4, 
F555W - F8UW < 0.4, and F606W - F8UW < 0.3), 
total stellar mass, and additional full-field dA v applied. 
In general the number of MS stars may be taken as an 
indication of confidence in the answer, but older popula- 
tions will also generally have fewer MS stars. As a result, 
this mapping is very approximate. The total stellar mass 
listed is highly dependent on the precise value selected 
for the IMF, and as such should not be taken as an exact 
measurement. It is intended only to compare relative am- 
plitudes of star formation between various SNR regions. 
We observe six fields that have no significant recent SF, 
and identify these as likely either Type la SNR or SNR 
resulting from runaway progenitor stars. They are not 
included in analysis of the distribution of recovered pro- 
genitor masses. 

To examine any underlying biases in our data selection, 
we plot the recovered ages against the 50% completeness 
magnitudes in Fig. [71 The top panel displays the bluer 
filter for each particular CMD, while the bottom is for 
the redder (F814W in all cases). Following application 
of our depth requirement, we find no evidence of any 
correlation between recovered age and depth of data. 

The top panel of Fig. [8] plots the recovered ages as a 
function of the number of MS stars detected. The bottom 
panel plots recovered ages as a function of dA v used. 
Note that these values only include full-field dA v , not 
the 0.5 mag of dA v included for all young populations. 
No correlation is observed in either comparison, which 
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suggests that our treatment of differential reddening from 
§2.4 does not bias us towards a specific age. 

While an in-depth discussion of all SNR is prohibitively 
long, it is worth exploring a representative sample of the 
various cases observed. In general, fields with fewer main 
sequence stars will tend to have larger uncertainties, and 
we list our results in Table [5] in order of the number 
of these stars. However, this generalization is not abso- 
lute, and this ordering should be taken as approximate. 
Finally, some of the regions display no young star forma- 
tion, which we list in a separate category. For illustrative 
purposes, we subjectively identify four classes of results. 
Below we examine four representative SNR for each case. 

3.1. K376: Obvious Young Star Formation 

In the left panel of Fig. |9] we display the color- 
magnitude diagram for the region around K376, plotted 
in red. We also plot our best-fit model CMD in the back- 
ground in grey scale, with darker regions corresponding to 
a larger expected stellar density. The data we used for 
K376 are from HST project 12055 (Brick 9, Field 14), 
and included ACS images in F475W and F814W. The 
photometry is fairly deep, with 50% completeness limits 
at F475W=27.1 and F814W=25.8. We find 6874 stars in 
the 50 pc region around the SNR, 350 of which are MS 
stars, making K376 a very well populated CMD. 

A qualitative inspection of the CMD indicates a bright, 
blue upper main sequence, indicative of a significant 
young stellar population. The CMD displays a great deal 
of differential reddening, as indicated by a large spread 
in the red clump. We used a full-field value of dA v = 0.7, 
in addition to dA v — 0.5 for young stars, to model this 
CMD. MATCH predictably finds that this CMD displays 
significant young star formation: we find the field well fit 
by a single burst of star formation at 12±1 Myr, corre- 
sponding to a mass estimate of 16±1 M The older burst 
at ^40 Myr is of significantly lower prominence compared 
to the young burst. The result is well-constrained due to 
the large number of stars in the field. 

Similar SNR include 2-020, 2-024, K934, K947, BW- 
60, BW-69, BW-74. 

3.2. K180: Well Defined Older Populations 

Fig. [TO] displays the CMD and SFH for the region 
surrounding K180. The data for K180 are from HST 
project 12073 (Brick 2, Field 11), and includes ACS im- 
ages in F475W and F814W. The photometry is quite 
deep, with 50% completeness limits of F475W=27.2 and 
F814W=25.9 respectively. We find 3950 stars in the 
CMD, 162 of which are MS stars. We used full-field 
value of dA v = 0.6, in addition to dA v = 0.5 for young 
stars, to model this CMD. 

The upper main sequence is not nearly as bright or 
prominent as that of K376, and as such we would pre- 
dict an older population with a lower mass progenitor 
star. MATCH agrees, finding the population is fit best 
by a single burst of star formation at 33±2 Myr, corre- 
sponding to a mass of 8.8±0.2 M©. As above, because of 
the large number of stars in the field, the mass is fairly 
well constrained. 

Similar SNR include 1-006, 1-008, 1-009, 1-010, 2-025, 
2-044, 2-046, 2-048, 2-050, K516, K526A, K574, K594, 
BW-18, BW-20, BW-31, BW-32, BW-76, BW-81, BW- 
84, BW-86, BW-102, and BW-110. 



3.3. K891: SF Spread Over a Wide Range 

Fig, im displays the CMD and SFH for the region sur- 
rounding K891. The data for K891 are from HST project 
12055, including ACS images in F475W and F814W. The 
photometry is again quite deep, with 50% completeness 
limits of F475W=27.6 and F814W=26.6 respectively. 
The CMD has 1441 total stars, including 149 MS stars. 
We used a full-field value of dA v = 1.1, in addition to 
dA v = 0.5 for young stars, to model this CMD. 

The CMD for K891 doesn't appear to comprise a single 
uniform population, and indeed has stars spread out at a 
wide variety of ages. MATCH finds two distinct star for- 
mation bursts: one in our youngest age bin at less than 
4.4 Myr, and another in an older bin at ^32 Myr. The 
error bars also indicate that the relative prominence of 
these bursts is such that we may not favor one burst over 
the other. While we find a median mass of ^9 M , the 
uncertainties allow this result to range over all almost 
masses past the minimum mass for core-collapse. K891 
is an example of a SNR where we have reasonable confi- 
dence in our answer, but this answer allows no constraint 
on the parameters we ultimately wish to measure. 

Similar SNR include 2-049, K446, K497, K525A, 
K527A, K856A, K908, K956A, BW-11, BW-39, BW-44, 
BW-61, BW-65, BW-66, BW-71, BW-77, BW-82, BW- 
89, BW-105, and BW-106. 

3.4. 2-028: No Recent Star Formation 

Fig. [12] displays the CMD and SFH for the region 
surrounding 2-028. The data for 2-028 is also from 
HST project 10273, and includes ACS images in F555W 
and F814W. The photometry has more shallow 50% 
limits than the CMDs above, with F555W=25.7 and 
F814W=25.4 respectively. The CMD has only 374 stars, 
a smaller number of detections than typical, including 
53 MS stars. We used a full-field value of dA v = 0.5, 
in addition to dA v = 0.5 for young stars, to model this 
CMD. 

The CMD clearly has a very sparse and dim main se- 
quence, leading us to predict that MATCH would find 
little to no young star formation. Indeed, MATCH finds 
the youngest star formation as occurring at older than 
90 Myr, well beyond the age of stellar populations that 
we would expect to produce progenitors of core-collapse 
SNe. We classify this SNR as the result of either a Type 
la SN or a runaway star that has left the young stellar 
population at its birth site. We do not use any results 
in this category as for purposes of analyzing the mass 
distribution. 

Similar SNR include 2-016, 2-021, 2-026, BW-19, and 
BW-36. 

4. DISCUSSION AND SUMMARY 

4.1. Indications of a Minimum Mass 

Our method measures the prominence of a burst as a 
fraction of the total recent SF. This approach requires 
that we define the period of time that we consider as 
"recent" SF. Specifically, we must identify the age range 
where star formation can produce stars massive enough 
to result in CCSNe. Theoretical arguments and obser- 
vational evidence point to a minimum mass necessary 
for progenitors to undergo core-collapse to a neutron 
star. Stars below this mass are generally assumed to 
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leave behind white dwarf stars, producing no SN ex- 
plosion. Thus one may constrain this cross-over mass 
by either measuring the maximum mass from which a 
star may create a white dwarf, or the minimum mass 
necessary for a star to explode. Measurements of white 
dwarfs have de fined a lower limit on this minimum mass 
of 6.3-7.1 M^ ([Williams et all [20093) . corresponding to 
an age of between ~55 and 63 Myr, and direct progeni- 
tor mass measurements has have converged on a va lue of 
8±1 M ? (|Smartt et alJl2009HBotticella et al.H2012D . cor- 
responding to an age of between ~33 and 55 Myr. 

In our observed progenitor mass distribution, we would 
expect to see the following behavior: above the mini- 
mum mass for core-collapse, we expect the distribution 
of progenitor masses to follow the IMF, assuming that 
the recent SFR is approximately constant. Below the 
minimum mass, the inferred progenitor mass should have 
no physical connection to the CCSNe process and should 
reflect random sampling of the SFR at >50 Myr, produc- 
ing an essentially flat distribution in inferred progenitor 
mass. 

In the left panel of Fig. Q21 we plot the distribu- 
tion of progenitor masses for a variety of assumed 
minimum masses. We vary the minimum mass from 
9.6 M©to 6.0 M ^Cour chosen m asses are mapped using 
the M zams from iMarigo et all (|2008ft and iGirardi et all 
(|2010f) isochrones). We find that for assumed minimum 
mass greater than 8.1 M Q , the distribution increases un- 
til the assumed minimum mass is reached. For assumed 
minimum mass values below 8.1 Mq, there is a peak 
between 7.5 and 8.5 Mq, below which the number of 
progenitors drops, suggesting a minimum mass in this 
range. 

To find the actual minimum mass, we lower the as- 
sumed minimum mass until the measured minimum mass 
no longer reflects this assumed value. We note the am- 
plitude of the peak in the distribution is greatest for an 
assumed minimum of 7.3 Mq. To quantify the location 
of this peak, we calculate the derivative of the number 
of progenitors as a function of mass. We assume that 
the maximum value of this derivative occurs at the min- 
imum mass a star undergoes core-collapse, as this value 
identifies the beginning of the peak. In the right panel 
of Fig. [T3]we plot the location of this maximum against 
our assumed minimum mass. We note that above an as- 
sumed minimum mass of 7.3 Mq, the peak value traces 
the assumed minimum mass. For 7.3 Mq and below, 
however, the peak value is always around ^7.5 Mq. If 
we assume our uncertainties are at least the width of the 
mass bins, the we find a minimum mass for core-collapse 
between 7.0 and 7.8 M Q . This ran ge is consistent with 
the observational m e asurem ents of ISmartt et al.1 ()2009l ) 
and iBotticella et al.1 ([20121 ) . We therefore have adopted 
the 44.7 to 50 Myr (7.7 to 7.3 M ) bin as the oldest 
included for our final results. 

4.2. Progenitor Mass Distribution 

In the left panel of Fig. [I4]we plot the histogram of me- 
dian progenitor masses, restricted to 7 Mq and above. In 
the right panel of Fig. [I4]we plot the cumulative fraction 
of progenitor masses. Unless otherwise noted, we assume 
a maximum mass of 120 M Q , although the choice of this 
value at high masses is essentially irrelevant to the overall 
distributions given the rarity of extremely massive stars. 



Qualitatively, the observed distribution shows a lack 
of the most massive stars when compared to a Salpctcr 
IMF (dN/dM oc M a , where a = -2.35). We per- 
formed a Kolomogorov-Smirnov (KS) test, assuming a 
single power law distribution. We found values of a out- 
side the range —2.7 > a > —4.4 inconsistent with the 
measured distribution at 95% confidence. Alternatively, 
we may consider a model distribution that is a Salpeter 
IMF (a = —2.35) up to some maximum mass, which we 
may vary. We found that this model was inconsistent 
with the data at 95% confidence at assumed maximum 
masses > 26 M Q . However, precise determination of this 
value is difficult in our survey due to the intrinsic rarity 
of massive stars and the mass spacing in our isochrincs 
(see Fig. [6]). Rather, this value represents the sort of 
mass range in which one must consider CCSN possible 
in order to maintain a Salpeter IMF. 

In either scenario, the full distribution of measured 
masses suggests that some fraction of massive stars are 
not exploding as CCSN. This result has interesting im- 
plications for CCSN physics. A wide variety of SN chan- 
nels have been explored both theoretically and observa- 
tionally in the literature. Theoretical predictions have 
explored the possibility of direct black hole formation 
beyond a certain mass threshold somewhere around ^25 
Mq . The manifestation of such events in an overall mass 
distribution would be an observed lack of progenitors 
beyond the mass threshold, essentially a more bottom- 
heavy IMF than that of all massive stars. While the 
reality is likely something more complicated than a well- 
defined threshold between CCSN and black hole forma- 
tion (many different scenarios likely combine to produce 
a complicated mass distribution), the qualitative effect 
w ill be that wh i ch we observe in our distribution. 

ISmartt et al.l <|2009f ) first identified the red supergiant 
problem, an observed lack of Type IIP progenitors be- 
tween 16 and 30 M Q . Many so lutions have been pro- 
posed to explain the problem (see iWalmswell fc Eldrided 
(|2012l) and references therein). In additio n, the recent 
SN 2012aw may fall in thi s mass range (jFraser et al.1 
120121: IVan Dvk et~aT1 l2012af ). suggesting the possibility 
that Type IIP progenitors do exist in this range and have 
simply not yet been observed in sufficient number. We 
identify six progenitors with median progenitor masses 
between 16 and 30 Mq, although the uncertainties on 
many of these are large. If our progenitors were sam- 
pled uniformly from a Salpeter IMF, we would expect 
to find ~10 progenitors (20%) in the mass range from 
16 to 30 Mq. Thus while we don't observe a complete 
lack of progenitors in the specified mass range, we do ob- 
serve fewer than we would expect given a Salpctcr IMF 
distribution. 

Finally, while we assumed our SNR catalogs consti- 
tuted a complete, unbiased sample, the possibility exists 
that selection effects in the catalogs lead to the lower 
end of the progenitor mass distribution being sampled 
more heavily. In particular, extremely massive progeni- 
tors are likely associated with strong H II regions, where 
identification of SNRs is a more difficult observational 
task. It is possible SNRs of this type are systematically 
undcrsampled in the survey. 

We attempted to quantify both the IMF slope and min- 
imum mass using Markov Chain Monte Carlo methods, 
but found that the data could not produce meaningful 
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constraints beyond those determined by our more sim- 
plistic techniques. We believe the chief reason for this 
is the size of errors due to differential reddening of the 
fields. Part of the problem is treatment of differential 
reddening as a top-hat distribution. In addition, the in- 
clination of M31 contributes to these high differential 
extinction values. The application of this technique to a 
less inclined galaxy would be of benefit in this respect. 
Finally, the simple addition of more progenitor mass es- 
timates would allow us to better constrain our analysis. 
We are currently performing identical analysis on an ad- 
ditional ^65 SNR in M33 in pursuit of these final two 
points. 

4.3. Type la Contamination 

Type la SNRs coincident with star-forming regions 
could in principle affect our mass distribution. To test 
the possibility of additional Type la contamination be- 
yond the 11% observed, we examined Galcx FUV fluxes 
at the sites of all SNR in our sample. We assumed that 
sites with the lowest FUV flux corresponded to possible 
older Type la sites which happened to be coincident with 
a small amount of recent SF. The eight lowest flux SNR 
(not including the 6 with no recent SF) were BW-18, 
BW-69, 1-006, 1-010, 2-024, 2-050, K891, and K956a. 
From Table [5J these additional SNRs in general have 
fewer MS stars and less total SF then most SNRs in the 
sample. We found that after removing these SNR, our 
observed distribution now ranged from —2.6 > a > —4.3, 
which is consistent with our earlier measurement. This 
suggests that additional Type la contamination has lit- 
tle effect on our overall result, and still results in an IMF 
that is steeper than Salpeter (-2.35). 

4.4. Summary 

Using resolved HST photometry, we have analyzed the 
stars surrounding 59 SNR in M31. Using CMD fitting, 
we calculate a SFH within a 50 pc radius of each SNR. 
We find that 53 of the SNR regions display significant ev- 
idence of recent star formation, which we use to age-date 
the progenitor star. The remaining six regions display no 
recent star formation, and we consider them either possi- 
ble Type la candidates or the result of massive runaway 
progenitor stars. 

We examine the distribution of progenitor masses for 
our CCSN candidates and find a lack of massive stars 
compared to a standard Salpeter IMF (dN/dM oc M a , 
where a = —2.35). If a uniform single IMF is as- 
sumed, we find values for a outside the range —2.7 > 
a > —4.4 inconsistent with the measured distribution 
at 95% confidence. Alternatively, if we consider a dis- 
tribution that is a Salpeter IMF up to some maximum 
mass, we place an upper limit on the maximum mass 
allowed at Muax ~ 26 Mq. We also estimate a mini- 
mum mass for core collapse of between 7.0 and 7.8 M Q , 
which is both gr eater than the maximu m mass for white 
dwarf collapse ([Williams et al.l l2009cf) and consistent 
with direct progenitor measurements (jSmartt "et~aT][200l 
iBotticella et alj|2012j ). 
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Figure 1. SN progenitor locations plotted on a star-subtracted Ha map of M31. SNRs in red correspond to those with no young SF 
detected. We consider these as possible Type la locations, or as the results of runaway stars. We note that a large amount of our CCSN 
candidates fall along the star forming arms of M31. This is clearly expected (CCSN candidates will tend to be found in regions of recent 
SF), but it also indicates the constraints imposed on our sample by the locations of archival HST data (star-forming regions are targeted 
significantly more than other areas). 
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Figure 2. Marigo ct al. ( 2008); Gira rdi et at] H2010I ) isochrones plotted in F555W against F555W - F814W for logarithmic age bins of 6.60 
through 8.00. Redder isochrones correspond to older ages. We adopt the 50 Myr isochrone as the maximum age for a CCSN progenitor 
star, corresponding to a Mzams °f 7.3 Mq. For clarity, we only plot every-other isochrone; our actual fitting procedure uses twice as many 
age bins over the same range. 
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Figure 3. Left Panel: Observed CMD for the region around BW-102, F814W plotted against F475W-F814W. In the background we plot 
the best-fit model CMD generated by MATCH in greyscale, assuming no additional differential reddening. The model expects a much 
brighter MS than observed at F475W < 22.5 and poorly models the stars at F475W ~ 25.5, FA75W - F814W ~ 1.75 when differential 
reddening is not accounted for. Right panel: Cumulative recent star formation history calculated for the corresponding CMD. This SFH 
was from a blind run using MATCH'S default values, with 0.5 magnitudes of differential reddening applied only to the youngest ages. Based 
on this SFH, we may only definite an upper limit on the age of 4.5 Myr or less, corresponding to a mass of greater than 52 Mq. The CMD 
displays no bright, blue upper MS stars, causing us to be suspicious of the extremely young progenitor star. In addition, the extended red 
clump indicates that differential reddening is not unique to the young stellar population. 



Final Isochrone Mass (M @ ) 

672619 15 13 11 10 9 8 




Age (Myr) 

Figure 4. Cumulative star formation history for BW-102. Different lines correspond to different values of field-wide differential reddening, 
with bluer lines corresponding to dA v =0.0 and redder lines corresponding to gL4„=1.5. This is in addition to the 0.5 magnitudes of dA v 
applied by default to the young stellar population. Note that some lines are not visible because they find the exact same age as other lines, 
and therefore show up underneath other results. The red line was found to have the best fit while still maintaining a value above that of 
foreground reddening. It corresponds to an age of 36 Myr and a mass of 8.5 Mq. 
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Figure 5. Left Panel: Wc plot the observed CMD of the region surrounding BW-102, plotted as magnitude in F475W vs. (F475W- 
F814W). We also plot the best-fit model created by MATCH in greyscale behind the observed CMD. The scale is printed to the right of 
the CMD. This model is found using the differential reddening procedure described in §2.4, leading us to use a value of dA v = 0.5. Right 
Panel: Cumulative star formation over the most recent 50 Myr. We plot the best fit (purple line) with high and low errors as calculated 
by the Monte Carlo tests (red lines). The cross-hatched highlighted region is corresponds to our 16% to 84% confidence interval, which we 
define in §2.3. We find the data best fit by a burst of SF at 36^^ Myr, corresponding to 8.5^q'j Mq (highlighted region). 




Figure 6. Marigo ct al. (2008); Girardi ct al. (2010) final isochronc masses plotted against isochrone age for metallicities of Z=0.004, 
Z=0.008 Z=0.019 (solar), and Z=0.030. Based on the known gas-phase metallicity of M31 and our own best-fit metallicity values, we adopt 
solar metallicities for all age-to-mass conversions in this paper. 
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Figure 7. Top Panel: Wc plot the recovered ages for the progenitor stars against the 50% completeness limit for the redder broadband 
filter used. Bottom Panel: Same as top, but for the bluer filter used. The filter selection is somewhat arbitrary, since the bluer filter may be 
F475W, F555W, or F606W depending on the dataset. The redder filter will always be F814W. The data points do not include results for 
which no young SF was found. Following our depth cut, we find no significant correlation between depth in either filter and the measured 
age. 
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Figure 8. Top panel: We plot the recovered ages for the progenitor stars against the number of MS stars detected in the region around the 
SNR. Bottom panel: We plot the recovered ages against the value of dA v used for ACS data only. These values include only the additional 
user-added dA v , not the 0.5 mag of dA v automatically included for young populations. The data points do not include results for which 
no young SF was found. We find no significant correlation between age and either dA v or number of stars. 
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Figure 9. Left Panel: In red we plot the observed CMD of the region surrounding K376, plotted as magnitude in F814W vs. (F475W- 
F814W). We also plot the best-fit model created by MATCH in greyscale behind the observed CMD. The scale is printed to the right 
of the CMD. A prominent upper main sequence is clearly visible in the CMD, indicating a young stellar population and a more massive 
progenitor. Right Panel: Cumulative star formation over the most recent 50 Myr. We plot the best fit (purple line) with high and low 
errors as calculated by the Monte Carlo tests (red lines). The cross-hatched highlighted region is corresponds to our 16% to 84% confidence 
interval, which we define in §2.5. We find the data best fit by a single burst of SF at 12±1 Myr, corresponding to 16±1 M© (highlighted 
region). 




F475W-F814W Age (Myr) 

Figure 10. Same as Fig. [9] but for K180. Left Panel: In red we plot the observed CMD of the region surrounding K180, plotted as 
magnitude in F814W vs. (F475W-F814W). We also plot the best-fit model created by MATCH in greyscale behind the observed CMD. 
The scale is printed to the right of the CMD. The main sequence is dimmer than that for K376, indicating an older population and thus 
a less massive progenitor. Right Panel: Cumulative star formation over the most recent 50 Myr. We plot the best fit (purple line) with 
high and low errors as calculated by the Monte Carlo tests (red lines). The cross-hatched highlighted region is corresponds to our 16% to 
84% confidence interval, which we define in §2.5. We find the data best fit by a single burst of SF at 33±2 Myr, corresponding to 8.8±0.2 
Mq (highlighted region). 
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Figure 11. Same as Fig. [9] but for K891. Left Panel: In red we plot the observed CMD of the region surrounding K891, plotted as 
magnitude in F814W vs. (F475W-F814W). We also plot the best-fit model created by MATCH in greyscale behind the observed CMD. The 
scale is printed to the right of the CMD. Right Panel: Cumulative star formation over the most recent 50 Myr. We plot the best fit (purple 
line) with high and low errors as calculated by the Monte Carlo tests (red lines). The cross-hatched highlighted region is corresponds to our 
16% to 84% confidence interval, which we define in §2.5. Unlike the previous examples, K891 seems to display two star formation events 
of similar prominence. We are unable to associate the progenitor star with either burst, and as such the only constraint we may offer is to 
say that the progenitor is younger than 36 Myr, corresponding to a mass of > 9 Mq. 




F555W-F814W Age (Myr) 

Figure 12. Same as Fig. [9] but for 2-028, with F555W instead of F475W. Left Panel: In red we plot the observed CMD of the region 
surrounding 2-028, plotted as magnitude in F814W vs. (F555W-F814W). We also plot the best-fit model created by MATCH in greyscale 
behind the observed CMD. The scale is printed to the right of the CMD. The main sequence is much dimmer than that for K376, indicating 
an older population and thus a less massive progenitor. Right Panel: Cumulative star formation over the most recent 50 Myr. We find no 
SF within the past 50 Myr for 2-028. We classify 2-028 as either a possible Type la remnant, or as a possible runaway star. 
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Figure 13. Left Panel: Plot of distribution of progenitor masses assuming different minimum mass for a star to undergo core collapse. 
We plot distributions for 6.0, 6.3, 6.6, 7.0, 7.3, 7.7, 8.1, 8.6, 9.0, and 9.6 Mq. We note the existence of a large peak of masses located at 
~8 Mq , which we argue represents the minimum core-collapse mass. We also note that this peak grows in size as we lower the minimum 
mass down to 7.3 Mq, where it is constant. This suggests that we include the full number of progenitors when we include all star formation 
down to 7.3 Mq, equivalent to the past 50 Myr of star formation. Right Panel: On the horizontal axis, we plot the minimum mass that 
we consider for purposes of determining recent star formation. On the vertical axis, we plot the corresponding mass for the largest value 
in the derivative of the progenitor mass distribution with respect to mass. We argue this spike represents the minimum mass for a star 
to undergo core-collapse. We find that above 7.3 Mq, this value follows the assumed minimum mass. For 7.3 Mgand below, this value 
remains at ~7.5 Mq. We assume that the precision with which we may define a value is ± one age bin, and that the error bars on this 
minimum mass correspond to the final isochrone mass for these bins. Using this analysis, we find a range of values for the minimum mass 
from 7.0 to 7.8 Mq. 
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Figure 14. Left Panel: Histogram of median progenitor masses below 52 Mq. Right Panel: Cumulative fraction of progenitor mass 
distribution. We overplot a reference Salpeter IMF. Using a KS-test, we find the cumulative mass distribution consistent with a power-law 
IMF of the form dN/dM oc M a with -2.7 > a > -4.4. We plot these two slopes, as well as a Salpeter IMF (dN/dM oc M -2 *). While 
the distribution of masses greater than 60 Mq is not shown, the fraction greater than 60 Mq is given by the value for the cumulative 
fraction at 60 Mq. 



20 



Table 1 List of SNR with Deep 2-Filter HST Data 



SNR ID 


RA (degrees) 


DEC (degrees) 


HST Field Project ID Instrument 


Filters w/ 50% Completeness Lii 








Maenier et al. 


(1995) SNR 








1-006 


10.6318 


41.1005 


POS-33 


10273 


WFPC-2 


F555W=24.5, 


F814W=23.4 


1-008 


10.7675 


41.6031 


POS-21 


10273 


ACS 


F555W=26.6, 


F814W=26.0 


1-009 


10.7975 


41.6256 


POS-23 


10273 


ACS 


F555W=25.6, 


F814W=25.3 


1-010 


10.7979 


41.4853 


POS-18 


10273 


WFPC-2 


F555W=24.9, 


F814W=23.7 


2-016 


10.3196 


40.9554 


G-87 


6671 


WFPC-2 


F555W=25.2, 


F814W=23.9 


2-020 


10.4508 


41.1138 


G-104 


10260 


ACS 


F606W=25.6, 


F814W=23.5 


2-021 


10.4773 


40.7866 


G-119 


6671 


WFPC-2 


F555W=26.2, 


F814W=25.0 


2-024 


10.5958 


41.0036 


POS-29 


10273 


ACS 


F555W=25.8, 


F814W=25.1 


2-025 


10.6448 


40.9688 


POS-29 


10273 


ACS 


F555W=25.9, 


F814W=25.1 


2-026 


10.6698 


41.0447 


POS-30 


10273 


ACS 


F555W=26.0, 


F814W=25.2 


2-028 


10.7376 


40.9698 


POS-41 


10273 


ACS 


F555W=25.7, 


F814W=25.4 


2-044 


11.1989 


41.4654 


B08-F10 


12075 


ACS 


F475W=27.3, 


F814W=26.0 


2-046 


11.2928 


41.5993 


B12-F17 


12071 


ACS 


F475W=27.5, 


F814W=26.1 


2-048 


11.3094 


41.6033 


B12-F17 


12071 


ACS 


F475W=27.5, 


F814W=26.1 


2-049 


11.3254 


41.8683 


B15-F08 


12056 


ACS 


F475W=27.4, 


F814W=26.0 


2-050 


11.3662 


41.8698 


B15-F07 


12056 


ACS 


F475W=27.4, 


F814W=26.1 








Williams ct al. 


(1995) SNR 








BW-11 


10.2467 


40.6081 


G-76 


11081 


WFPC-2 


F606W=26.2, 


F814W=24.9 


BW-18 


10.3988 


41.1155 


G-104 


10260 


ACS 


F606W=26.1, 


F814W=24.9 


BW-19 


10.5408 


40.9472 


POS-27 


10273 


ACS 


F555W=26.1, 


F814W=25.5 


BW-20 


10.5433 


40.8644 


POS-26 


10273 


ACS 


F555W=26.0, 


F814W=25.0 


BW-31 


10.7317 


40.9956 


POS-41 


10273 


ACS 


F555W=26.4, 


F814W=25.4 


BW-32 


10.7329 


40.9717 


POS-41 


10273 


ACS 


F555W=25.7, 


F814W=24.6 


BW-36 


10.7725 


41.3750 


G-205 


10260 


ACS 


F606W=25.5, 


F814W=24.5 


BW-39 


10.7933 


41.6282 


POS-23 


10273 


ACS 


F555W=25.8, 


F814W=25.4 


BW-44 


10.8779 


41.6882 


POS-24 


10273 


ACS 


F555W=26.5, 


F814W=25.0 


BW-60 


11.0883 


41.9018 


B15-F12 


12056 


ACS 


F475W=28.0, 


F814W=26.8 


BW-61 


11.1054 


41.3501 


B06-F10 


12105 


ACS 


F475W=27.3, 


F814W=25.9 


BW-65 


11.1488 


41.4227 


B06-F04 


12105 


ACS 


F475W=27.3, 


F814W=25.9 


BW-66 


11.1550 


41.8666 


B15-F17 


12056 


ACS 


F475W=27.5, 


F814W=26.2 


BW-69 


11.1825 


41.9645 


B17-F18 


12059 


ACS 


F475W=27.6, 


F814W=26.5 


BW-71 


11.1958 


41.4886 


B08-F04 


12075 


ACS 


F475W=27.5, 


F814W=26.0 


BW-74 


11.2129 


41.4847 


B08-F04 


12075 


ACS 


F475W=27.1, 


F814W=25.8 


BW-76 


11.2267 


41.5121 


B08-F04 


12075 


ACS 


F475W=27.3, 


F814W=25.9 


BW-77 


11.2269 


41.5306 


B08-F04 


12075 


ACS 


F475W=27.3, 


F814W=25.9 


BW-81 


11.2858 


41.6101 


B12-F17 


12071 


ACS 


F475W=27.4, 


F814W=26.1 


BW-82 


11.2892 


41.8523 


B15-F08 


12056 


ACS 


F475W=27.6, 


F814W=26.3 


BW-84 


11.3162 


41.6561 


B12-F11 


12071 


ACS 


F475W=27.4, 


F814W=26.0 


BW-86 


11.3387 


41.6668 


B12-F11 


12071 


ACS 


F475W=27.3, 


F814W=26.0 


BW-89 


11.3650 


41.9036 


B15-F01 


12056 


ACS 


F475W=27.5, 


F814W=26.3 


BW-102 


11.4675 


42.1618 


B21-F11 


12055 


ACS 


F475W=27.6, 


F814W=26.8 


BW-105 


11.6296 


41.9886 


B18-F03 


12108 


ACS 


F475W=27.6, 


F814W=26.6 


BW-106 


11.6417 


42.1804 


B21-F08 


12055 


ACS 


F475W=28.0, 


F814W=27.0 


BW-110 


11.6896 


42.2183 


B21-F01 


12055 


ACS 


F475W=27.7, 


F814W=27.0 



Braun fc Walterbos (1993) SNR 



K180 


10.9186 


41.1814 


B02-F11 


12073 


ACS 


F475W= 


=27.2, 


F814W= 


=25.9 


K376 


11.0850 


41.5804 


B09-F14 


12057 


ACS 


F475W= 


=27.1, 


F814W= 


=25.8 


K446 


11.1298 


41.3572 


B06-F04 


12105 


ACS 


F475W= 


=27.3, 


F814W= 


=26.0 


K497 


11.1562 


41.4133 


B06-F10 


12105 


ACS 


F475W= 


=27.3, 


F814W= 


=26.0 


K516/BW-67 


11.1699 


41.4145 


B06-F04 


12105 


ACS 


F475W= 


=27.4, 


F814W= 


=26.0 


K525A 


11.1820 


41.4372 


B08-F10 


12075 


ACS 


F475W= 


=27.3, 


F814W= 


=26.0 


K526A 


11.1715 


41.4653 


B08-F10 


12075 


ACS 


F475W= 


=27.3, 


F814W= 


=26.0 


K527A 


11.1834 


41.4465 


B08-F10 


12075 


ACS 


F475W= 


=27.3, 


F814W= 


=26.0 


K574 


11.2101 


41.4649 


B08-F10 


12075 


ACS 


F475W= 


=27.4, 


F814W= 


=26.1 
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Table 1 List of SNR with Deep 2-Filter HST Data 



SNR ID 


RA (degrees) 


DEC (degrees) 


HST Field 


Project ID 


Instrument 


Filters w/ 50% Completeness Limits 


K594 


11.2201 


41.9161 


B15-F10 


12056 


ACS 


F475W=27.2, F814W=26.0 


K856A 


11.4313 


41.9313 


B16-F05 


12106 


ACS 


F475W=27.6, F814W=26.4 


K891 


11.5405 


42.2198 


B21-F04 


12055 


ACS 


F475W=27.6, F814W=26.6 


K908 


11.6231 


41.9685 


B18-F03 


12108 


ACS 


F475W=28.0, F814W=27.0 


K934/BW-107 


11.6467 


42.2266 


B21-F01 


12055 


ACS 


F475W=27.6, F814W=26.9 


K947/2-047 


11.6689 


42.1911 


B21-F07 


12055 


ACS 


F475W=27.6, F814W=26.8 


K956A 


11.6792 


42.2171 


B21-F01 


12055 


ACS 


F475W=27.7, F814W=27.0 



Table 2 Progenitor Age and Mass Results 



SNR ID 


A/zams(M ) 


Age (Myr) 


# MS Stars 


Total Stars 


Total Mass Formed (1O 2 M ) 


Additional dA 


BW-74 


n+f 

— z 


11+? 


825 


4572 


169 




0.3 


BW-44 


> 10 


< 26 


6 11 


2085 


3 




1.5 


BW-86 


9 6+ 29 


™±li 


575 


4344 


IS 




0.3 


BW-84 
K527A 


VTtf 


34±? 2 
lit? 


572 
531 


5158 
4925 


19 
62 




0.3 
0.1 


2-049 
BW-81 


76 +34 

' -°-0.3 

l2+\ 


19t| 


517 
502 


4364 
5212 


60 
16 




0.6 
0.4 


K525A 
BW-65 


20±? 2 
« 1+ 7 - 6 

°- 1 -0.40 


n 4+35 

40+^ 


501 
500 


4677 
5263 


121 
57 




0.2 

0.5 


BW-77 




47±| 7 


473 


4968 


31 




0.2 


K934 


>13 


>16 


455 


1934 


75 




0.7 


K908 
BW-66 


8-613:2 

36l 6 27 


5 5+ 26 
°- -0.5 


418 
411 


3360 
3310 


17 
63 




0.3 
1.6 


BW-71 
BW-31 
K446 


lOtf 
9 9+ - 3 
9.0tg 6 4 


25+ 19 
— 20 

27+. 2 

32+ 4 
26 


401 
397 
389 


5281 
2288 
4921 


37 
23 
26 




0.6 
0.3 
0.6 


K594 


8.4+ ,; 2 


38+J 


379 


2537 


26 




0.5 


BW-106 


19±?i 


lOlf 


372 


2642 


61 




1.3 


K497 
K856A 


10±J 


7+43 
' -1 


360 
356 


5106 
2125 


26 
54 




0.2 

0.3 


K376 


i6±i 


12±J 


350 


6874 


55 




0.7 


BW-76 


o c+2.5 
8 - b -0.3 


37t? 4 


345 


5013 


35 




0.3 


1-008 
BW-61 


15±f 


46l 18 
13±J 8 


336 
303 


1895 
4835 


61 
106 




1.7 
0.8 


2-050 
K526A 
BW-82 
BW-89 

K947 
BW-60 


7- 6to.3 
8 4+ - 2 

o i+8 

8 - 1 -0A 

i3if 


47 -24 

38± 2 
40t* 7 

40l| 6 

16+ 2 

I5±l 


289 
288 
281 
265 
245 
243 


4343 
4842 
4017 
4296 
1874 
3277 


11 

3 
8 
13 
37 
39 




0.5 
0.3 
0.3 
0.9 
0.9 
1.5 


BW-105 


>8 


<44 


239 


1945 


13 




0.3 


BW-11 

1- 009 
BW-20 

2- 044 
2-046 


8 9+ - 8 
7 9+0.3 

nti 


S3+. 3 , 

42t 3 3 
o 7 +i3 
a ' -12 

2ll' 


235 
219 
192 
186 
185 


1199 
739 
1339 
4368 
4445 


105 
36 
22 
13 
16 




0.9 
1.2 
1.8 
0.2 
0.1 


2-024 


16±i 


i2±i 


164 


1934 


5 




0.5 


BW-39 


>11 


<20 


163 


844 


104 




2.1 


K574 
K180 
2-048 
K891 


o i+6 
8a -0.4 

« H+0- 2 
8.8_ 2 

— U.4 

>9 


40ll 4 
33+j 
45+1, 
<36 


162 
162 
160 
149 


4207 
3950 
3581 
1441 


18 
23 
33 
23 




0.5 
0.6 
0.7 
1.1 


2-025 
BW-18 


9 9+ 05 


39+i 1 
27+j 


140 
117 


1287 
2165 


13 
25 




0.3 
2.0 


K516 
K956A 
BW-69 


8 - 5 io.3 

8-lto 8 4 
™±% 


36+ 4 
00 -14 

40^33 

7.5+_l 


109 
103 
100 


2210 
1457 
1467 


10 
22 

8 




0.4 
1.1 
1.5 


BW-32 
BW-110 


9 3+ - 3 
12±l 


30± 2 
20±1 


91 

86 


562 
1391 


6 
12 




0.5 
0.8 



Supernova Remnant Progenitor Masses in M31 



Table 2 Progenitor Age and Mass Results 



SNR ID 




Age (Myr) # MS Stars Total Stars 


Total Mass Formed(10 2 M Q ) 


Additional dA v Applied 


BW-102 


^"^—0.4 


00 11 1 


1568 


31 


1.5 


2-020 


18+1 


ll j OO 


1137 


20 


1. 1 


1-006 


L1 -l 


24+2 10 


601 


6 


0.0 


1-U1U 


19 +1 
12 -1 


19+5 10 


335 


9 


n 




Fields with no 


recent SF (Probable Type la 


or Runaway Progenitors), not included in 


distribution 


BW-19 




291 


2520 





0.5 


BW-36 




10 


4162 





1.3 


2-026 




191 


2421 





0.1 


2-021 




122 


1268 





0.0 


2-028 




53 


374 





0.5 


2-016 




5 


466 





0.6 



